<html>
 <head>
  <meta charset="utf-8"/>
  <meta content="width=device-width, initial-scale=1, maximum-scale=1, user-scalable=no" name="viewport"/>
  <title>
   ZhangDixiang: R软件在精算教学中的应用案例  | 数螺 | NAUT IDEA
  </title>
  <link href="http://cdn.bootcss.com/bootstrap/3.3.6/css/bootstrap-theme.min.css" rel="stylesheet"/>
  <link href="http://cdn.bootcss.com/bootstrap/3.3.6/css/bootstrap.min.css" rel="stylesheet"/>
  <style type="text/css">
   #xmain img {
                  max-width: 100%;
                  display: block;
                  margin-top: 10px;
                  margin-bottom: 10px;
                }

                #xmain p {
                    line-height:150%;
                    font-size: 16px;
                    margin-top: 20px;
                }

                #xmain h2 {
                    font-size: 24px;
                }

                #xmain h3 {
                    font-size: 20px;
                }

                #xmain h4 {
                    font-size: 18px;
                }


                .header {
	           background-color: #0099ff;
	           color: #ffffff;
	           margin-bottom: 20px;
	        }

	        .header p {
                  margin: 0px;
                  padding: 10px 0;
                  display: inline-block;  
                  vertical-align: middle;
                  font-size: 16px;
               }

               .header a {
                 color: white;
               }

              .header img {
                 height: 25px;
              }
  </style>
  <script src="http://cdn.bootcss.com/jquery/3.0.0/jquery.min.js">
  </script>
  <script src="http://cdn.mathjax.org/mathjax/latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML" type="text/javascript">
   MathJax.Hub.Config({
          tex2jax: {inlineMath: [['$','$'], ['\\(','\\)']]},
          TeX: {equationNumbers: {autoNumber: "AMS"}}
        });
  </script>
  <script src="http://nautstatic-10007657.file.myqcloud.com/static/css/readability.min.js" type="text/javascript">
  </script>
  <script type="text/javascript">
   $(document).ready(function() {
                 var loc = document.location;
                 var uri = {
                  spec: "http://cos.name/2011/01/the-application-of-r-in-actual-science-with-case-study/",
                  host: "http://cos.name",
                  prePath: "http://cos.name",
                  scheme: "http",
                  pathBase: "http://cos.name/"
                 };
    
                 var documentClone = document.cloneNode(true);
                 var article = new Readability(uri, documentClone).parse();
     
                 document.getElementById("xmain").innerHTML = article.content;
                });
  </script>
  <!-- 1466441558: Accept with keywords: (title(0.25):案例,门户网站,教学,统计学,服务平台,中国,软件,精算, topn(0.4):概率,费用率,男表,统计学,R语言,数值,教学,文件,缔香,均值,精算,参数,计算,命名,死亡率,保险,向量,赔款,函数,状态,生命表,矩阵,命令,寿险,拟合,年龄,案例,数据,费率,软件).-->
 </head>
 <body class="single single-post postid-2869 single-format-standard sidebar" onload="">
  <div class="header">
   <div class="container">
    <div class="row">
     <div class="col-xs-6 col-sm-6 text-left">
      <a href="/databee">
       <img src="http://nautidea-10007657.cos.myqcloud.com/logo_white.png"/>
      </a>
      <a href="/databee">
       <p>
        数螺
       </p>
      </a>
     </div>
     <div class="hidden-xs col-sm-6 text-right">
      <p>
       致力于数据科学的推广和知识传播
      </p>
     </div>
    </div>
   </div>
  </div>
  <div class="container text-center">
   <h1>
    ZhangDixiang: R软件在精算教学中的应用案例
   </h1>
  </div>
  <div class="container" id="xmain">
   <div class="hfeed site" id="page">
    <header class="site-header" id="masthead" role="banner">
     <div id="cos-logo">
      <a href="http://cos.name/">
       <img src="http://cos.name/wp-content/themes/COS-kermesinus/images/headers/cos-logo.png"/>
      </a>
     </div>
     <div class="navbar" id="navbar">
      <nav class="navigation main-navigation" id="site-navigation" role="navigation">
       <h3 class="menu-toggle">
        菜单
       </h3>
       <a class="screen-reader-text skip-link" href="http://cos.name/2011/01/the-application-of-r-in-actual-science-with-case-study/#content" title="跳至内容">
        跳至内容
       </a>
       <div class="menu-%e6%88%91%e7%9a%84%e8%8f%9c%e5%8d%95-container">
        <ul class="nav-menu" id="menu-%e6%88%91%e7%9a%84%e8%8f%9c%e5%8d%95">
         <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-home menu-item-4746" id="menu-item-4746">
          <a href="http://cos.name">
           主页
          </a>
         </li>
         <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-has-children menu-item-8120" id="menu-item-8120">
          <a href="http://cos.name/cn">
           论坛
          </a>
          <ul class="sub-menu">
           <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-8124" id="menu-item-8124">
            <a href="http://cos.name/cn/wp-login.php?action=register">
             论坛注册
            </a>
           </li>
           <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-8125" id="menu-item-8125">
            <a href="http://cos.name/cn/wp-login.php">
             论坛登录
            </a>
           </li>
          </ul>
         </li>
         <li class="menu-item menu-item-type-post_type menu-item-object-page menu-item-8110" id="menu-item-8110">
          <a href="http://cos.name/books/">
           图书资料
          </a>
         </li>
         <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-8112" id="menu-item-8112">
          <a href="http://cos.name/videos">
           视频教程
          </a>
         </li>
         <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-8333" id="menu-item-8333">
          <a href="http://cos.name/salon/">
           统计沙龙
          </a>
         </li>
         <li class="menu-item menu-item-type-post_type menu-item-object-page menu-item-8111" id="menu-item-8111">
          <a href="http://cos.name/chinar/">
           R语言会议
          </a>
         </li>
         <li class="menu-item menu-item-type-post_type menu-item-object-page menu-item-8109" id="menu-item-8109">
          <a href="http://cos.name/training/">
           讲座与培训
          </a>
         </li>
         <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-8115" id="menu-item-8115">
          <a href="http://cos.name/cn/forum/comprehensive/job/">
           招聘信息
          </a>
         </li>
         <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-has-children menu-item-4780" id="menu-item-4780">
          <a href="http://cos.name/about">
           关于我们
          </a>
          <ul class="sub-menu">
           <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-8113" id="menu-item-8113">
            <a href="http://cos.name/2008/11/how-to-work-with-cos/">
             加入我们
            </a>
           </li>
           <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-8119" id="menu-item-8119">
            <a href="http://cos.name/donate/">
             赞助我们
            </a>
           </li>
           <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-8114" id="menu-item-8114">
            <a href="http://cos.name/cn/forum/26">
             项目合作
            </a>
           </li>
          </ul>
         </li>
        </ul>
       </div>
       <form action="http://cos.name/" class="search-form" method="get" role="search">
        <label>
         <span class="screen-reader-text">
          搜索：
         </span>
        </label>
       </form>
      </nav>
      <!-- #site-navigation -->
     </div>
     <!-- #navbar -->
    </header>
    <!-- #masthead -->
    <div class="site-main" id="main">
     <div class="content-area" id="primary">
      <div class="site-content" id="content" role="main">
       <article class="post-2869 post type-post status-publish format-standard hentry category-actuary tag-r-language tag-20 tag-300" id="post-2869">
        <header class="entry-header">
         <h1 class="entry-title">
          R软件在精算教学中的应用案例
         </h1>
         <div class="entry-meta">
          <span class="date">
           <a href="http://cos.name/2011/01/the-application-of-r-in-actual-science-with-case-study/" rel="bookmark" title="链向R软件在精算教学中的应用案例的固定链接">
            <time class="entry-date" datetime="2011-01-12T23:22:04+00:00">
             2011/01/12
            </time>
           </a>
          </span>
          <span class="categories-links">
           <a href="http://cos.name/category/data/actuary/" rel="category tag">
            风险精算
           </a>
          </span>
          <span class="tags-links">
           <a href="http://cos.name/tag/r-language/" rel="tag">
            R语言
           </a>
           、
           <a href="http://cos.name/tag/%e5%ad%a6%e4%b9%a0%e7%bb%8f%e5%8e%86/" rel="tag">
            学习经历
           </a>
           、
           <a href="http://cos.name/tag/%e7%b2%be%e7%ae%97/" rel="tag">
            精算
           </a>
          </span>
          <span class="author vcard">
           <a class="url fn n" href="http://cos.name/author/fatsmellycat/" rel="author" title="查看所有由张缔香发布的文章">
            张缔香
           </a>
          </span>
         </div>
         <!-- .entry-meta -->
        </header>
        <!-- .entry-header -->
        <div class="entry-content">
         <p>
          <strong>
           本文作者为张缔香，文章由COS编辑部审核发表，略有修改。
           <a href="http://cos.name/wp-content/uploads/2011/01/R软件在精算教学中的应用案例.pdf" target="_blank">
            点击此处下载/阅读本文PDF版本
           </a>
          </strong>
         </p>
         <hr/>
         R软件做为一种统计软件，因其开源、免费、灵活的诸多优点得到越来越多的关注，无论网络上还是实体书店，关于R的教程铺天盖地，不甚枚举。因此，本文的目标不是做R的教程，而是将R和保险、精算教学结合起来，通过几个案例来说明R在保险、精算专业日常的教学和研究中可用之处。
         <p>
         </p>
         <p>
          作者在保险、精算的理论、专业知识方面水平有限，也没有编程基础，使用R也只2个月左右的时间，所以以下几个案例，无论在理论还是R编程技巧方面都有所欠缺。但做为加入COS的首篇文章，当求真求实，不做文抄公。
         </p>
         <p>
          最后，本文的几个案例，均来自于作者在南开大学保险系读硕期间的一些课程作业，主要目的在于完成某一项任务，而没有过多地考虑通用性、用户界面友好性等，比如所有代码中几乎都没有写入报错、警告语句等。当然，作者也假定本文读者不单熟悉R的常用语法、函数等，而且具有相应的统计学、保险、精算专业知识，案例中的术语、结论等都直接引用，不做解释。
         </p>
         <h2 style="text-align: center;">
          案例一 简单寿险产品的保费、准备金等
         </h2>
         <p>
          《南开大学经济学系列实验教材》丛书里有一本李秀芳老师主编《寿险精算实务实验教程》，是李秀芳老师2010年秋学期主讲的《寿险精算实务实验》课程的指定教材，作者也选修了这一门课。因此本案例主要来自这门课的实验内容，但在实现方法上，不再使用换算函数法，而转向概率的思路，利用精算等价原理来求得寿险产品的净保费、毛保费、准备金等。
         </p>
         <p>
          本案例中的R脚本文件命名为li.R，存放生命表的CSV文件命名为clt03.csv。
         </p>
         <h3>
          1. 寿险计价函数 li()
         </h3>
         <p>
          <strong>
           描述：
          </strong>
          以列表的形式返回目标寿险产品的毛保费、理论责任准备金、费用。
         </p>
         <p>
          <strong>
           参数：
          </strong>
         </p>
         <ul>
          <li>
           投保人年龄age：单个数值型，取值范围0-105；
          </li>
          <li>
           生存给付bf.s：年初被保险人生存条件下的给付，非负数值型向量；
          </li>
          <li>
           死亡给付bf.d：被保险人死亡的年度末的给付，非负数值型向量；
          </li>
          <li>
           定价利率ir：按年计量的复利，数值型向量；
          </li>
          <li>
           费用率feerate：预定费用占毛保费的比例，非负数值型向量；
          </li>
          <li>
           保费费率prerate：均衡、递增、增减、或是其他形式的，数值型向量;
          </li>
          <li>
           定价生命表file：字符串型对象，存放生命表的.CSV文件的文件路径；
          </li>
          <li>
           产品类型type：取值范围{1、2、3、4}，分别对应男性、女养老金业务，男性、女性非养老金业务。
          </li>
         </ul>
         <p>
          <strong>
           注释：
          </strong>
         </p>
         <h5>
          1. 生存给付、死亡给付的设置方式
         </h5>
         <p>
          以保额1000元，20年定期死亡保险为例，那么死亡给付 bf.d=rep(1000,20)，生存给付 bf.s=0；
          <br/>
          若是保额1000元的20年定期年金，则 bf.d=0，bf.s=rep(1000,20)；
          <br/>
          若保险给付不是有规律的序列，则必须以 c(x,y,…) 的形式具体地给出。
         </p>
         <h5>
          2. 费用率、保费费率的设置方式
         </h5>
         <p>
          费用率、保费费率必须设置为数值型向量，且长度必须与目标寿险产品的缴费期限相同。
          <br/>
          以10年缴费为例，假定第一年的费用率0.6，第二年0.4，第三年0.2，第四年及以后为0.1，那么feerate=c(0.6,0.4,0.2,rep(0.1,7))。若费用率为0，则必须指明为0 ，不可省略。
          <br/>
          同样以10年缴费为例，若是均衡保费，即每一年的保费都相同，则保费费率 prerate=rep(1,10)。递增型 prerate=c(10:1) 表示第一年缴10单位保费，第二年9单位，以此类推。
         </p>
         <h5>
          3. 生命表 CSV 文件的格式
         </h5>
         <p>
          生命表部分数据如下所示，存放在 CSV 文件里：
         </p>
         <pre class="brush: plain">age	 male	    female      malep	     femalep
0	0.000722	0.000661	0.000627	0.000575
1	0.000603	0.000536	0.000525	0.000466
2	0.000499	0.000424	0.000434	0.000369
3	0.000416	0.000333	0.000362	0.00029
4	0.000358	0.000267	0.000311	0.000232
5	0.000323	0.000224	0.000281	0.000195</pre>
         <p>
          4. 本案例中采用的生命表的范围是0至105岁，所以在设置年龄、生存、死亡给付、费用率、保费费率时，年龄和其它四个参数长度最大者之和不能超过106。
          <br/>
          <strong>
           示例：
          </strong>
         </p>
         <p>
          目标产品：20岁投保、20年死亡保险、定价利率2.5%、10年均衡缴费、男性非养老金业务的净保费
         </p>
         <pre class="brush: r">source(“li.R”)
li(age=20,bf.s=0,bf.d=rep(1000,20),ir=0.025,feerate=rep(0,10),prerate=rep(1,10),file="clt03.csv",type=1)</pre>
         <p>
          输出结果如下：
         </p>
         <pre class="brush: r">$premium
 [1] 1.623953 1.623953 1.623953 1.623953 1.623953 1.623953 1.623953
 [8] 1.623953 1.623953 1.623953</pre>
         <pre class="brush: r">$reserve
 [1]  1.036481  2.059095  3.076657  4.096164  5.119730  6.148513
 [7]  7.183697  8.229508  9.282238 10.335088  9.710075  9.017654
[13]  8.244942  7.390757  6.447828  5.406646  4.256400  2.982906
[19]  1.569502  0.000000  0.000000  0.000000  0.000000  0.000000</pre>
         <pre class="brush: r">$fee
 [1] 0 0 0 0 0 0 0 0 0 0</pre>
         <p>
          <strong>
           附录：
          </strong>
          附录一
         </p>
         <h3>
          2. 应用案例：
         </h3>
         <p>
          做为《寿险实务实验课程》的结课作业，作者设计了一款少儿保险产品，并用R和li() 函数计算净保费、毛保费、现金价值、法定责任准备金、利润测试等，由于篇幅过大，待整理完毕之后再发上来。
         </p>
         <h2 style="text-align: center;">
          案例二 BMS奖罚系统
         </h2>
         <p>
          本案例的R脚本文件命名为bms.R，索赔次数的CSV格式数据文件命名为BMSclaim.csv。
         </p>
         <h3>
          1. 目标BM系统
         </h3>
         <p>
          本案例的目标BMS系统采用中国保险行业协会2007年4月制定的奖罚系统，具体如下所示：
         </p>
         <table align="center" frame="hsides">
          <caption>
           表1 中国保险行业协会制奖罚系统表(2007.04)
          </caption>
          <tbody>
           <tr>
            <th>
             规则
            </th>
            <th>
             保费等级系数
            </th>
            <th>
             等级
            </th>
           </tr>
           <tr>
            <th>
             连续3年及3年以上无赔款
            </th>
            <th>
             <span style="font-weight: normal;">
              0.7
             </span>
            </th>
            <th>
             <span style="font-weight: normal;">
              1
             </span>
            </th>
           </tr>
           <tr>
            <th>
             连续2年无赔款
            </th>
            <th>
             <span style="font-weight: normal;">
              0.8
             </span>
            </th>
            <th>
             <span style="font-weight: normal;">
              2
             </span>
            </th>
           </tr>
           <tr>
            <th>
             上年无赔款
            </th>
            <th>
             <span style="font-weight: normal;">
              0.9
             </span>
            </th>
            <th>
             <span style="font-weight: normal;">
              3
             </span>
            </th>
           </tr>
           <tr>
            <th>
             新保或上年发生3次以下赔款
            </th>
            <th>
             <span style="font-weight: normal;">
              1.0
             </span>
            </th>
            <th>
             <span style="font-weight: normal;">
              4
             </span>
            </th>
           </tr>
           <tr>
            <th>
             上年发生3次赔款
            </th>
            <th>
             <span style="font-weight: normal;">
              1.1
             </span>
            </th>
            <th>
             <span style="font-weight: normal;">
              5
             </span>
            </th>
           </tr>
           <tr>
            <th>
             上年发生4次赔款
            </th>
            <th>
             <span style="font-weight: normal;">
              1.2
             </span>
            </th>
            <th>
             <span style="font-weight: normal;">
              6
             </span>
            </th>
           </tr>
           <tr>
            <th>
             上年发生5次及5次以上赔款
            </th>
            <th>
             <span style="font-weight: normal;">
              1.3
             </span>
            </th>
            <th>
             <span style="font-weight: normal;">
              7
             </span>
            </th>
           </tr>
          </tbody>
         </table>
         <h3>
          2. 转移概率矩阵
         </h3>
         <p style="text-align: center;">
          $
          <br/>
          M=\left( \begin{array}{ccc} {1-\sum_{i=0}^4p_i} \; p_4 \; p_3 \; p_1+p_2 \; p_0 \; 0 \; 0 \\ {1-\sum_{i=0}^4p_i} \; p_4 \; p_3 \; p_1+p_2 \; p_0 \; 0 \; 0 \\ {1-\sum_{i=0}^4p_i} \; p_4 \; p_3 \; p_1+p_2 \; p_0 \; 0 \; 0 \\ {1-\sum_{i=0}^4p_i} \; p_4 \; p_3 \; p_1+p_2 \; p_0 \; 0 \; 0 \\ {1-\sum_{i=0}^4p_i} \; p_4 \; p_3 \; p_1+p_2 \; 0 \; p_0 \; 0 \\ {1-\sum_{i=0}^4p_i} \; p_4 \; p_3 \; p_1+p_2 \; 0 \; 0 \; p_0 \\ {1-\sum_{i=0}^4p_i} \; p_4 \; p_3 \; p_1+p_2 \; 0 \; 0 \; p_0\\ \end{array} \right)
          <br/>
          $
         </p>
         <p>
          其中$p_k=\frac{\lambda^k}{k!}e^{-\lambda}$，表示投保人在一个保险年度中恰好发生 k 次索赔的概率。
         </p>
         <h3>
          3. 计算稳态分布
         </h3>
         <p>
          在稳定状态的处于各个级别的概率$\pi=(\pi_6,\pi_7,\dots,\pi_1)$满足：
         </p>
         <p style="text-align: center;">
          $\sum_i\pi_i=1$, \ $\pi=\pi M$
         </p>
         <p>
          解得：
         </p>
         <p style="text-align: center;">
          $\pi_7=1-\sum_0^4p_i; \ \pi_6=p_4; \ \pi_5=p_3; \ \pi_4=p_2+p_1; \ \pi_3=p_0-p_0^2; \ \pi_2=p_0^2-p_0^3; \ \pi_1=p_0^3$
         </p>
         <h3>
          4. 利用索赔次数的实际数据求出具体的转移概率矩阵
         </h3>
         <p>
          本案例的索赔次数的观察值数据以如下所示，存放在CSV文件里。
         </p>
         <pre class="brush: plain">claim	exposure
0	   27141
1	   5789
2	   1443
3	   457
4	   155
5	   56
6	   27
7	   2
8	   2
9	   1
10	   0</pre>
         <p>
          这里采用泊松逆高斯分布来拟合索赔次数分布，采用极大似然估计的方法来估计参数(在做这个案例时，作者对R了解甚少，很多的功能、函数其实是可以在各种专业package里找到的，这里不详述。)。结果可通过以下的命令查看：
         </p>
         <pre class="brush:r">source("bms.R")
gamma_pig
beta_pig</pre>
         <p>
          按此估计的参数值，恰好发生k次索赔的概率可通过以下的命令查看：
         </p>
         <pre class="brush: r">print(pr)</pre>
         <p>
          按以上的索赔概率值，和事先设定好的转移规则，可得到转移概率矩阵，通过以下命令查看：
         </p>
         <pre class="brush: r">tmat</pre>
         <p>
          为求稳定状态下处于各保费等级的概率值，将转移矩阵取30次方(实际上只需要5次方左右即可达到稳定的极限分布)，这个极限矩阵的各行已经趋于相同，也即所求的稳定状态下的保费等级分布。可能通过以下命令查看极限分布：
         </p>
         <pre class="brush: r">tmat_lim</pre>
         <h3>
          5. 计算评价BMS系统优劣的指标
         </h3>
         <h5>
          5.1 稳定状态下的平均保费率（新保费率为1.0）
         </h5>
         <p>
          稳定状态下的平均保费率即是以稳定状态下处于各保费等级的概率分布为权重，各状态下的保费费率的加权平均值，计算结果为0.8196。可以通过以下命令查看：
         </p>
         <pre class="brush: r">avp</pre>
         <h5>
          5.2 相对稳定状态平均水平(RSAL)
         </h5>
         <p>
          RSAL=(稳定平均水平 – 最低水平)/(最高水平 – 最低水平)=0.1993，可以通过以下命令查看：
         </p>
         <pre class="brush: r">rsal</pre>
         <h5>
          5.3 对新司机的隐性惩罚(ECL)
         </h5>
         <p>
          ECL=(新保费率 – 稳定平均费率)/稳定平均费率=0.2201，可以通过以下命令查看：
         </p>
         <pre class="brush: r">ecl</pre>
         <h5>
          5.4 变异系数COD=标准差/均值
         </h5>
         <p>
          各年的平均保费以及变异系数可以通过以下命令查看：
         </p>
         <pre class="brush: r">cod_mat</pre>
         <p>
          平均保费和变异系数逐年变化曲线可以通过以下的命令查看：
         </p>
         <pre class="brush: r">graph1()</pre>
         <h5>
          5.5 风险区分度
         </h5>
         <p>
          由于对索赔频率的数据采用的是泊松逆高斯拟合的，所以每个保费等级内不同保单的索赔频率服从参数为$\lambda$的泊松分布(条件分布)；在所有保费等级中$\lambda$又服从逆高斯分布，其参数的极大似然估计已在前文交待。
         </p>
         <p>
          对于各保费等级的保单，其索赔频率和均值和方差即为以上的复合分布的无条件均值和方差。均值和方差的计算均要计算无穷积分。在R里，采用离散化的方式来近似计算，即采取足够大的区间、足够小的间隔，通过求和的方式来计算无穷积分的近似值。各保费等级的内的索赔频率和方差的估计值以及相应结构 的参数可以通过以下命令来查看：
         </p>
         <pre class="brush: r">fin_mat</pre>
         <p>
          各个保费等级下的索赔频率曲线图可以通过以下的命令查看：
         </p>
         <pre class="brush: r">graph2()</pre>
         <h5>
          5.6 BMS的弹性
         </h5>
         <p>
          BMS的弹性公式为：$elac=(dp/p)/(d^{\lambda}/\lambda)=d \ln p /d \ln \lambda$，该公式为连续形式下的弹性公式，在实验中稳定状态下平均保费的函数形式的精确表达式难以计算出，只能采取离散化的方式来近期计算。具体思路就是把微分用差分代替，取间隔够小、数量够多的$\lambda$及对应的稳定状态下的平均保费，得到弹性的曲线图。可以通过以下的命令查看：
         </p>
         <pre class="brush: r">graph3()</pre>
         <p>
          <strong>
           代码：
          </strong>
          附录二
         </p>
         <h2 style="text-align: center;">
          案例三 利用LeeCarter模型进行死亡率预测
         </h2>
         <p>
          本案例中的脚本文件命名为forecast.R，两个CSV文件分别命名为93pm.csv 和 03pm.csv。
         </p>
         <h3>
          1. 数据来源及初步整理
         </h3>
         <p>
          选用中国生命表90-93养老金男表和00-03养老金男表做为本案例的基础数据。把101至105岁划分为一组，100岁及其以上每一整数年龄分为一组，共102组。100岁及其以上组死亡率直接采用死亡率qx，100岁以上组采用中心死亡率。采用线性插值法得出1994-2002年间的死亡率。90-93养老金男表和00-03养老金男表各存放于一个CSV文件中，按如下图所示的格式存放数据：
         </p>
         <p>
          在R的命令窗口输入以下命令来观察1993-2003 死亡率矩阵：
         </p>
         <pre class="brush: r">source(“forecast.R”)
tab9303</pre>
         <p>
          把死亡率取自然对数，然后减去每一年龄的死亡率的对数在1993-2003年间的均值，得到残差矩阵A。在R的命令窗口输入以下命令来观察均值和A：
         </p>
         <pre class="brush: r">ax
amat</pre>
         <p>
          <a href="http://cos.name/wp-content/uploads/2011/01/actual-fig-1.png">
           <img src="http://cos.name/wp-content/uploads/2011/01/actual-fig-1-300x300.png"/>
          </a>
          <br/>
          各年龄的均值如上图所示。
         </p>
         <h3>
          2. SVD分解
         </h3>
         <h5>
          2.1 k值
         </h5>
         <p>
          利用 R 里的 svd() 命令对死亡率的残差矩阵A做SVD分解。通过以下命令可观察分解得出的U、S、V矩阵：
         </p>
         <pre class="brush: r">svda$u
svda$d
svda$v</pre>
         <p>
          把S矩阵右乘V的转置矩阵，取其结果的第一行，即为Lee-Carter 模型中的kt。按1993-2003的顺序，kt的各个值如下：
         </p>
         <p>
          -2.2891315， -1.9135574， -1.5181148， -1.1002924， -0.6570094，-0.1844099，
          <br/>
          0.3224553， 0.8701444， 1.4677250， 2.1284910， 2.8736996
         </p>
         <p>
          可以在R里输入以下命令观察kt的值及其总和：
         </p>
         <pre class="brush:r">kt
sum(kt)</pre>
         <h5>
          2.2 b值
         </h5>
         <p>
          在SVD分解得出的U矩阵里，取其第一列，即可得到未经中心化的bx的值。本人采用保险研究增刊里的方法，采用线性拟合的方法来求bx的值。在R里输入以下命令可观察bx的拟合值以及拟合结果：
         </p>
         <pre class="brush: r">bmat</pre>
         <p>
          在bmat矩阵中其中，如各列名称所示，第一列为年龄，第二列为bx，第三至六列分别为线性拟合bx时的t值、t临界值、F值、F临界值和拟合优度值，可以看出所有的bx都是显著的。
          <br/>
          <a href="http://cos.name/wp-content/uploads/2011/01/actual-fig-2.png">
           <img src="http://cos.name/wp-content/uploads/2011/01/actual-fig-2-300x300.png"/>
          </a>
         </p>
         <h5>
          2.3 k值的调整
         </h5>
         <p>
          采用最小化RMSPE的方法来得到最优的kt值，RMSPE的公式如下：
         </p>
         <p>
          采用非线性最小值的算法来求解使得RMSPE最小的kt值，具体如下：
         </p>
         <p>
          -2.3445627， -1.9458174， -1.5271026， -1.0869498， -0.6237423， -0.1356942
         </p>
         <p>
          0.3791744， 0.9230672， 1.4984474， 2.1080823， 2.7550978
         </p>
         <p>
          通过如下命令来观察求出的k值：
         </p>
         <pre class="brush: r">kt_adj</pre>
         <h5>
          2.4 预测
         </h5>
         <p>
          对kt值和时间变量1-11(1993年做为1，2003做为11，或者直接用1993，1994，……2003区别只在系数的值，对预测没有影响)，得出kt的向后15期的预测值如下：
         </p>
         <p>
          3.044373， 3.551768， 4.059164， 4.566559， 5.073955， 5.581350， 6.088746， 6.596141， 7.103537， 7.610932， 8.118328， 8.625723， 9.133118， 9.640514， 10.147909
         </p>
         <p>
          可以通过如下命令查看kt的15期预测值及其95%置信区间：
         </p>
         <pre class="brush:r">k.lm</pre>
         <p>
          使用以上的kt预测值来预测后15年的出生平均余命，
         </p>
         <p>
          80.95256 ，81.35479 ，81.74651 ，82.12812， 82.49997 ，82.86240， 83.21572，83.56023，
         </p>
         <p>
          83.89621， 84.22392 ，84.54362， 84.85554， 85.15991， 85.45694， 85.74685
         </p>
         <p>
          下图为向后预测39期得出的出生平均余命的预测值，在2050年左右，中国男性的平均出生余命达到90岁。这似乎有些过高，所以本模型只适合做15年以内的预测。
         </p>
         <p>
          <a href="http://cos.name/wp-content/uploads/2011/01/actual-fig-3.png">
           <img src="http://cos.name/wp-content/uploads/2011/01/actual-fig-3-300x300.png"/>
          </a>
         </p>
         <p>
          <strong>
           代码：
          </strong>
          附录三
         </p>
         <h3>
          附录
         </h3>
         <p>
          [box]
         </p>
         <p>
          <a href="http://cos.name/wp-content/uploads/2011/01/R软件在精算教学中的应用案例（附录）.zip">
           R软件在精算教学中的应用案例（附录）
          </a>
          （包含附录一、二、三）
         </p>
         <p>
          [/box]
         </p>
         <div class="wumii-hook">
          <br/>
          <br/>
         </div>
        </div>
        <!-- .entry-content -->
        <footer class="entry-meta">
         <div class="author-info">
          <div class="author-avatar">
           <img src="http://sdn.geekzu.org/avatar/1f9f8e6dd55cf77702f453eebb74e2dd?s=74&amp;d=monsterid&amp;r=g"/>
          </div>
          <!-- .author-avatar -->
          <div class="author-description">
           <h2 class="author-title">
            关于张缔香
           </h2>
           <p class="author-bio">
            Statistics@XJTU
Insurance&amp;Actuarial science@NKU
            <a class="author-link" href="http://cos.name/author/fatsmellycat/" rel="author">
             查看所有由张缔香发表的文章
             <span class="meta-nav">
              →
             </span>
            </a>
           </p>
          </div>
          <!-- .author-description -->
         </div>
         <!-- .author-info -->
        </footer>
        <!-- .entry-meta -->
       </article>
       <!-- #post -->
       <nav class="navigation post-navigation" role="navigation">
        <h1 class="screen-reader-text">
         文章导航
        </h1>
        <div class="nav-links">
         <a href="http://cos.name/2011/01/cache-objects-in-sweave-stat-computation-and-graphics/" rel="prev">
          <span class="meta-nav">
           ←
          </span>
          Sweave后传：统计报告中的大规模计算与缓存
         </a>
         <a href="http://cos.name/2011/01/publishing-promotion-and-collaboration-in-statistics/" rel="next">
          统计学论文的发表流程、及统计学家的晋升和合作（内幕）
          <span class="meta-nav">
           →
          </span>
         </a>
        </div>
        <!-- .nav-links -->
       </nav>
       <!-- .navigation -->
       <div class="comments-area" id="comments">
        <h2 class="comments-title">
         《
         <span>
          R软件在精算教学中的应用案例
         </span>
         》有1个想法
        </h2>
        <ol class="comment-list">
         <li class="comment even thread-even depth-1" id="comment-7055">
          <article class="comment-body" id="div-comment-7055">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/be7b469e7f4c8485a287e63ff490cb5f?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              tjdina
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2011/01/the-application-of-r-in-actual-science-with-case-study/#comment-7055">
              <time datetime="2016-01-28T12:35:27+00:00">
               2016/01/28 12:35
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             请问张老师，您的“R软件在精算教学中的应用案例（附录）（包含附录一、二、三）”及“代码：附录三”为什么不能下载？
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给tjdina" class="comment-reply-link" href="http://cos.name/2011/01/the-application-of-r-in-actual-science-with-case-study/?replytocom=7055#respond" onclick='return addComment.moveForm( "div-comment-7055", "7055", "respond", "2869" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
         </li>
         <!-- #comment-## -->
        </ol>
        <!-- .comment-list -->
        <div class="comment-respond" id="respond">
         <h3 class="comment-reply-title" id="reply-title">
          发表评论
          <small>
           <a href="http://cos.name/2011/01/the-application-of-r-in-actual-science-with-case-study/#respond" id="cancel-comment-reply-link" rel="nofollow" style="display:none;">
            取消回复
           </a>
          </small>
         </h3>
         <form action="http://cos.name/wp-comments-post.php" class="comment-form" id="commentform" method="post" novalidate="">
          <p class="comment-notes">
           <span id="email-notes">
            电子邮件地址不会被公开。
           </span>
           必填项已用
           <span class="required">
            *
           </span>
           标注
          </p>
          <p class="comment-form-comment">
           <label for="comment">
            评论
           </label>
           <textarea aria-required="true" cols="45" id="comment" maxlength="65525" name="comment" required="required" rows="8">
           </textarea>
          </p>
          <p class="comment-form-author">
           <label for="author">
            姓名
            <span class="required">
             *
            </span>
           </label>
          </p>
          <p class="comment-form-email">
           <label for="email">
            电子邮件
            <span class="required">
             *
            </span>
           </label>
          </p>
          <p class="comment-form-url">
           <label for="url">
            站点
           </label>
          </p>
          <p class="form-submit">
          </p>
          <p style="display: none;">
          </p>
          <p style="display: none;">
          </p>
         </form>
        </div>
        <!-- #respond -->
       </div>
       <!-- #comments -->
      </div>
      <!-- #content -->
     </div>
     <!-- #primary -->
     <div class="sidebar-container" id="tertiary" role="complementary">
      <div class="sidebar-inner">
       <div class="widget-area">
        <aside class="widget widget_text" id="text-3">
         <h3 class="widget-title">
          关注统计之都
         </h3>
         <div class="textwidget">
          <ul>
           <li>
            新浪微博
            <a href="http://weibo.com/cosname">
             @统计之都
            </a>
           </li>
           <li>
            人人网
            <a href="http://renren.com/cosname">
             @统计之都
            </a>
           </li>
           <li>
            Twitter
            <a href="http://twitter.com/cos_name">
             @cos_name
            </a>
           </li>
          </ul>
         </div>
        </aside>
        <aside class="widget widget_text" id="text-6">
         <h3 class="widget-title">
          微信公众平台
         </h3>
         <div class="textwidget">
          <img src="http://cos.name/wp-content/uploads/2013/04/qrcode-8cm.jpg"/>
          <p style="font-size:12px;margin-left:15px">
           微信号 CapStat
          </p>
          <p>
           我们将第一时间向您推送主站和论坛的精彩内容，以及统计之都的线下活动、竞赛、培训和会议信息。
          </p>
         </div>
        </aside>
        <aside class="widget widget_text" id="text-8">
         <h3 class="widget-title">
          站内导航
         </h3>
         <div class="textwidget">
          <ul>
           <li class="page_item">
            <a href="http://cos.name/cn/">
             中文论坛
            </a>
           </li>
           <li class="page_item">
            <a href="http://cos.name/videos/">
             视频
            </a>
           </li>
           <li class="page_item">
            <a href="http://cos.name/salon/">
             统计沙龙
            </a>
           </li>
           <li class="page_item">
            <a href="http://cos.name/books/">
             图书出版
            </a>
           </li>
           <li class="page_item">
            <a href="http://cos.name/training/">
             教育培训
            </a>
           </li>
           <li class="page_item">
            <a href="http://cos.name/about/">
             关于我们
            </a>
           </li>
           <li class="page_item">
            <a href="http://cos.name/donate/">
             捐赠
            </a>
           </li>
           <li class="page_item">
            <a href="http://cos.name/chinar/">
             R语言会议
            </a>
           </li>
          </ul>
         </div>
        </aside>
        <aside class="widget widget_links" id="linkcat-2">
         <h3 class="widget-title">
          友情链接
         </h3>
         <ul class="xoxo blogroll">
          <li>
           <a href="http://stat.ruc.edu.cn" target="_blank" title="中国人民大学统计学院网站">
            中国人民大学统计学院
           </a>
          </li>
          <li>
           <a href="http://rucdmc.net">
            中国人民大学数据挖掘中心
           </a>
          </li>
          <li>
           <a href="http://birc.gsm.pku.edu.cn/" target="_blank">
            北京大学商务智能研究中心
           </a>
          </li>
          <li>
           <a href="http://sam.cufe.edu.cn/" target="_blank" title="中央财经大学统计与数学学院网站">
            中央财经大学统计与数学学院
           </a>
          </li>
          <li>
           <a href="http://tjx.cueb.edu.cn/" target="_blank" title="首都经济贸易大学统计学院网站">
            首经贸统计学院
           </a>
          </li>
          <li>
           <a href="http://www.shookr.com/">
            数客网大数据社区
           </a>
          </li>
          <li>
           <a href="http://www.xueqing.tv/" target="_blank" title="数据科学在线学习平台">
            雪晴数据网
           </a>
          </li>
          <li>
           <a href="http://iera.name/" target="_blank" title="IERA是一个旨在普及、传播和增进工业工程知识的非营利性网站">
            IERA（直通IE）
           </a>
          </li>
         </ul>
        </aside>
        <aside class="widget widget_categories" id="categories-2">
         <h3 class="widget-title">
          全部分类
         </h3>
         <label class="screen-reader-text" for="cat">
          全部分类
         </label>
         <select class="postform" id="cat" name="cat">
          <option value="-1">
           选择分类目录
          </option>
          <option class="level-0" value="925">
           cos访谈  (4)
          </option>
          <option class="level-0" value="659">
           图书出版  (5)
          </option>
          <option class="level-0" value="379">
           数学方法  (14)
          </option>
          <option class="level-1" value="381">
           分析与代数  (1)
          </option>
          <option class="level-1" value="380">
           概率论  (9)
          </option>
          <option class="level-1" value="382">
           随机过程  (5)
          </option>
          <option class="level-0" value="210">
           数据分析  (81)
          </option>
          <option class="level-1" value="203">
           多元统计  (3)
          </option>
          <option class="level-1" value="42">
           数据挖掘与机器学习  (42)
          </option>
          <option class="level-1" value="36">
           生物与医学统计  (17)
          </option>
          <option class="level-1" value="35">
           计量经济学  (4)
          </option>
          <option class="level-1" value="296">
           金融统计  (3)
          </option>
          <option class="level-1" value="303">
           风险精算  (7)
          </option>
          <option class="level-0" value="177">
           模型专题  (15)
          </option>
          <option class="level-1" value="38">
           回归分析  (10)
          </option>
          <option class="level-1" value="41">
           时间序列  (2)
          </option>
          <option class="level-0" value="784">
           每周精选  (24)
          </option>
          <option class="level-1" value="183">
           可视化  (9)
          </option>
          <option class="level-0" value="967">
           沙龙纪要  (3)
          </option>
          <option class="level-0" value="18">
           经典理论  (46)
          </option>
          <option class="level-1" value="37">
           抽样调查  (3)
          </option>
          <option class="level-1" value="4">
           统计推断  (26)
          </option>
          <option class="level-1" value="236">
           试验设计  (7)
          </option>
          <option class="level-1" value="39">
           非参数统计  (3)
          </option>
          <option class="level-0" value="1">
           统计之都  (279)
          </option>
          <option class="level-1" value="884">
           中国R会议  (2)
          </option>
          <option class="level-1" value="885">
           中国R语言会议  (5)
          </option>
          <option class="level-1" value="446">
           出国留学  (3)
          </option>
          <option class="level-1" value="179">
           推荐文章  (90)
          </option>
          <option class="level-1" value="3">
           新闻通知  (75)
          </option>
          <option class="level-1" value="263">
           统计刊物  (10)
          </option>
          <option class="level-1" value="174">
           网站导读  (40)
          </option>
          <option class="level-1" value="204">
           职业事业  (51)
          </option>
          <option class="level-1" value="213">
           高校课堂  (9)
          </option>
          <option class="level-0" value="178">
           统计计算  (28)
          </option>
          <option class="level-1" value="40">
           优化与模拟  (15)
          </option>
          <option class="level-1" value="43">
           贝叶斯方法  (6)
          </option>
          <option class="level-0" value="378">
           软件应用  (116)
          </option>
          <option class="level-1" value="44">
           统计图形  (36)
          </option>
          <option class="level-1" value="110">
           统计软件  (83)
          </option>
         </select>
        </aside>
        <aside class="widget widget_recent_comments" id="recent-comments-2">
         <h3 class="widget-title">
          最新评论
         </h3>
         <ul id="recentcomments">
          <li class="recentcomments">
           <span class="comment-author-link">
            fineboom
           </span>
           发表在《
           <a href="http://cos.name/2016/06/use-shiny-fleetly-set-up-visual-prototype-system/#comment-7317">
            利用shiny包快速搭建可视化原型系统
           </a>
           》
          </li>
          <li class="recentcomments">
           <span class="comment-author-link">
            胡家新
           </span>
           发表在《
           <a href="http://cos.name/2016/06/r%e8%af%ad%e5%8d%83%e5%af%bb%e7%ac%ac%e4%b8%89%e6%9c%9f%ef%bc%9a%e5%bc%a0%e6%97%a0%e5%bf%8c%e7%a9%b6%e7%ab%9f%e7%88%b1%e8%b0%81%ef%bc%9f/#comment-7316">
            R语千寻第三期：张无忌究竟爱谁？
           </a>
           》
          </li>
          <li class="recentcomments">
           <span class="comment-author-link">
            <a class="url" href="http://www.zijiacha.com/category.php?id=6" rel="external nofollow">
             南糯山普洱茶
            </a>
           </span>
           发表在《
           <a href="http://cos.name/2016/06/r%e8%af%ad%e5%8d%83%e5%af%bb%e7%ac%ac%e4%b8%89%e6%9c%9f%ef%bc%9a%e5%bc%a0%e6%97%a0%e5%bf%8c%e7%a9%b6%e7%ab%9f%e7%88%b1%e8%b0%81%ef%bc%9f/#comment-7315">
            R语千寻第三期：张无忌究竟爱谁？
           </a>
           》
          </li>
          <li class="recentcomments">
           <span class="comment-author-link">
            J
           </span>
           发表在《
           <a href="http://cos.name/2016/05/value-of-the-reputation-from-the-data/#comment-7314">
            数据告诉你：高信誉的卖家应该收高价，还是收低价？
           </a>
           》
          </li>
          <li class="recentcomments">
           <span class="comment-author-link">
            <a class="url" href="http://gg" rel="external nofollow">
             Hilda
            </a>
           </span>
           发表在《
           <a href="http://cos.name/2013/01/drawing-map-in-r-era/#comment-7311">
            R时代，你要怎样画地图？
           </a>
           》
          </li>
         </ul>
        </aside>
        <aside class="widget widget_rss" id="rss-282869971">
         <h3 class="widget-title">
          <a class="rsswidget" href="http://cos.name/cn/topics/feed/">
           <img src="http://cos.name/wp-includes/images/rss.png"/>
          </a>
          <a class="rsswidget" href="http://cos.name/cn/topics/feed/">
           中文论坛新帖
          </a>
         </h3>
         <ul>
          <li>
           <a class="rsswidget" href="http://cos.name/cn/topic/417372/">
            处理时间数据和产生时间序列的问题
           </a>
          </li>
          <li>
           <a class="rsswidget" href="http://cos.name/cn/topic/6790/">
            《统计陷阱》下载 （How to lie with statistics）
           </a>
          </li>
          <li>
           <a class="rsswidget" href="http://cos.name/cn/topic/818/">
            统计学的世界（第五版）
           </a>
          </li>
          <li>
           <a class="rsswidget" href="http://cos.name/cn/topic/16574/">
            class(x) 返回值值是AsIs,AsIs代表什么，有什么用处？
           </a>
          </li>
          <li>
           <a class="rsswidget" href="http://cos.name/cn/topic/417366/">
            如何用R绘制一个分类算法的决策规则
           </a>
          </li>
         </ul>
        </aside>
        <aside class="widget widget_text" id="text-2">
         <h3 class="widget-title">
          登录/RSS
         </h3>
         <div class="textwidget">
          <ul>
           <li>
            <a href="http://cos.name/wp-admin/">
             登录
            </a>
           </li>
           <li>
            <a href="http://cos.name/feed/" title="使用 RSS 2.0 同步站点内容">
             文章
             <abbr title="Really Simple Syndication">
              RSS
             </abbr>
            </a>
           </li>
           <li>
            <a href="http://cos.name/comments/feed/" title="RSS 上的最近评论">
             评论
             <abbr title="Really Simple Syndication">
              RSS
             </abbr>
            </a>
           </li>
          </ul>
         </div>
        </aside>
       </div>
       <!-- .widget-area -->
      </div>
      <!-- .sidebar-inner -->
     </div>
     <!-- #tertiary -->
    </div>
    <!-- #main -->
    <footer class="site-footer" id="colophon" role="contentinfo">
     <div class="site-info">
      版权所有 © 2014 统计之都 | 由
      <a href="http://wordpress.org/">
       WordPress
      </a>
      构建 | 主题修改自
      <a href="http://wordpress.org/themes/twentythirteen">
       Twenty Thirteen
      </a>
     </div>
     <!-- .site-info -->
    </footer>
    <!-- #colophon -->
   </div>
   <!-- #page -->
   <p style="margin:0;padding:0;height:1px;overflow:hidden;">
    <a href="http://www.wumii.com/widget/relatedItems" style="border:0;">
     <img src="http://static.wumii.cn/images/pixel.png"/>
    </a>
   </p>
  </div>
 </body>
</html>